Wave packet dynamics in hole Luttinger systems 
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For hole systems with an effective spin 3/2 we analyzed analytically and numerically the evolution 
of wave packets with the different initial polarizations. The dynamics of such systems is determined 
by the 4x4 Luttinger Hamiltonian. We work in the space of arbitrary superposition of light- and 
heavy-hole states of the "one-particle system". For 2D packets we obtained the analytical solution 
for the components of wave function and analyzed the space-time dependence of probability densities 
as well as angular momentum densities. Depending on the value of the parameter a — kod (fco is the 
average momentum vector and d is the packet width) two scenarios of evolution are realized. For 
a >> 1 the initial wave packet splits into two parts and the coordinates of packet center experience 
the transient oscillations or Zitterbeuiegung (ZB) as for other two-band systems. In the case when 
a « 1 the distribution of probability density at t > remains almost cylindrically symmetric and 
the ripples arise at the circumference of wave packet. The ZB in this case is absent. We evaluated 
and visualized for different values of parameter a the space-time dependence of angular momentum 
densities, which have the multipole structure. It was shown that the average momentum components 
can precess in the absence of external or effective magnetic fields due to the interference of the light- 
and heavy hole states. For localized initial states this precession has a transient character. 

PACS numbers: 73.22.-f, 73.63.Fg, 78.67.Ch, 03.65.Pm 



I. INTRODUCTION 

In a crystalline solids described by the Hamiltonian 

H = p /2m + V(r), where p is the momentum operator, 
V{r) is a periodic potential, and m is the electron mass in 
vacuum, the electron velocity operator v = p/m does not 
commute with the Hamiltonian H and so, the velocity is 
not a constant of motion. As a result, various approxima- 
tions including nearly-free, tightly bound electron mod- 
els or k ■ p methods lead to unusual dynamics of charged 
particles, which is completely different from their average 
behavior. For the motion of electrons in solids, a highly 
oscillatory component appears in the time evolution of 
physical observables such as position, velocity and spin 
angular momentum. This phenomenon which is known 
as Zitterbeuiegung (ZB) was described by Schrodinger in 
1930 for a free relativistic electron in vacuumi^ It was 
later understood that the ZB (a trembling motion) is due 
to the interference of states with positive and negative 
electron energies. 

At first the ZB phenomena in crystalline solids was 
predicted with the use of LCAO method in Refs.[3-5]. 
Later an oscillatory motion of electron wave packets has 
been considered in wide class of 3D solids and nanos- 
tructures, including narrow gap semiconductors 6 , car- 
bon nanotubes^, 2D electron gas with Rashba spin-orbit 
coupling 8,9 , single and bilayer graphene 10 ' 11 , and also 
superconductors^. The ZB of photons near the Dirac 
point in 2D photonic crystals was discussed in Ref[13]. 
Rusin and Zavadski considered ZB in crystalline solids 
for nearly-free and tightly bound electrons and concluded 
that ZB in solids is a rule rather than an exception.— 
They determined the parameters of trembling motion 
and concluded that, when the bands are decoupled, elec- 



trons should be treated as particles of a finite size. In 
another work by Rusin and Zavadski was performed the 
calculations of ZB in the Kronig-Penney model and was 
demonstrated that the two-band (k-p) model is adequate 
for this description^ 

The authors of Ref. [16] studied the semiclassical 
motion of holes by numerical solution of equations for 
Heisenberg operators in the Luttinger model at the pres- 
ence of dc electric field. The trajectories and angular 
momentum oscillations of heavy and light holes reminis- 
cent of the ZB were found and analyzed. 

In Refs. [17,18] was presented theoretical analysis of 
ZB for systems characterized by different Hamiltonians 
with gapped or spin-orbit coupled spectrum, including 
Luttinger, Kane, Rashba and Dresselhaus. It was demon- 
strated the analogy of ZB in all these systems and pre- 
sented a unified treatment of these phenomena. Culccr 
et al. showed that in the Luttinger model the time de- 
pendence of the angular momentum operator S(t) corre- 
sponds to a spin precession in the absence of any exter- 
nal or effective magnetic fields In Ref. [20] the formal- 
ism for treating the wave packet evolution in solids with 
two band energy spectrum has been developed, where 
the nontrivial non-Abelian terms in the equation of mo- 
tion, arising from the additional degree of freedom, were 
obtained. In particular, the effect of spin separation in 
electric field was analyzed. 

Mainly the ZB like effects were considered in the frame 
of Heisenberg representation. In this representation the 
time evolution of the system can be inferred directly from 
the differential equations that governed the behavior of 

operators r(t), v(t) and S(t). As usual these studies have 
dealt with the behavior of expectation values of these 
operators without going into the details of the space-time 
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evolution. 

In this work we study the space-time evolution of the 
wave packets in Schrodinger representation for the Lut- 
tinger model which describes the states of heavy and light 
holes in vicinity of the top of energy band in a wide class 
of p - doped III-V semiconductors. We obtain not only 
the average values such as average velocities or average 
angular momentum but more informative characteristics 
of hole wave packets: the probability density and effective 
spin densities. We demonstrate that space-time dynamics 
of wave packets has, as a rule, a complex character and 
it helps us to explain many interesting patterns of evolu- 
tion for different initial parameters and polarizations of 
wave packets. 

The paper is organized as follows. In section II the ana- 
lytical expressions for the 4-component wave function are 
obtained for different initial angular momentum polariza- 
tions. Two different scenarios of the space-time evolution 
of probability densities for different values of the param- 
eter a — kod, where fco is average momentum and d is the 
initial packet width are discussed. The transient oscilla- 
tions of the position operator are analyzed. The effective 
spin dynamics is considered in Sec. III. We evaluate and 

visualized the effective spin densities and average S(t). 
It was shown that the average angular momentum ex- 
periences transient precession about the average packet 
momentum kg when S(0) is not perpendicular or parallel 
to fco- We conclude with some remarks in Sec. IV. 
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The two-fold degenerate energy eigenvalues of H for 
light holes (L) and heavy holes {H ) are 



£L(fc)~(-7i + 2 72 )fc 2 , 

h 2 

E H {k) = - — {~ n + 2 l2 )k 2 . 
Am 



(2) 



The helicity operator defined by A = commutes 
with Hamiltonian (1) so that a helicity A (A = ±^,±|) 
is a good quantum number: A = ±^ corresponds to the 
light holes and A = ±| to the heavy holes. 

In this work we consider 2D model with fc = (k x , k y ) 
that allows us to obtain the simple analytical results 
concerning the time-evolution of the initial wave packet 
^>(x,y). In this case the eigenstates of the helicity oper- 
ator ^\ are 



II. SPACE - TIME EVOLUTION AND 

ZITTERBE WEGUNG 



Holes near the top of valence band of common semi- 
conductors such as Ge and GaAs have an effective spin 
j = 3/2 and are described by the 4x4 Luttinger 
Hamiltonian^, which in the isotropic approximation has 



A=±5 




H = 



2m 



(-7i + -y 2 )k 2 - 2 l2 (Sk) 2 



(1) 



Here 71 and 72 are Luttenger parameters, m is the elec- 
tron mass in vacuum, and S is the vector of 4 x 4 spin 
matrices correspondent to spin S = |. The matrices S x , 
S y , S z in Eq.(l) are given by 



S x — i 
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Here we will use more simple eigenfunctions of the 
Hamiltonian (1) which are the linear combinations of 



^L(H),i{r) = n^f) U L{H),i, « = 1,2, (3) 
where ifi^(r) = e lkr /27r and spinors U^^h)^ are given by 



Uli 



HI 




, U H 2 




(4) 
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with k x — k cos ip, k y — kshiip. 

Now the most general wave function can be written as 



^{r,t)= f dk n (f)lc L 4k)U L:l e 



iE L t/h, 



+CH,i{k)U H>i e 



-iE H t/H 



(5) 



where C^rm i(k) are to be determined from the Fourier 
expansion of &(r) at t = 0. 

i) As a first example we compute the time evolution of 
the Gaussian spinor 



(6) 




FIG. 1: (Color online) The full hole density for initial wave 
packet, Eq.(6) at the moment t = 0.2 (in the units of to = 
for 71 = 7.65 and 72 = 2.41, for width d — 3 ■ 10 _6 cra 
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and fc = 2 ■ 10 6 cm 



It realizes a wave packet with positive average momen- 
tum and the effective spin oriented initially along z axis, 
S z = |. The appropriate coefficients can be readily found 
from Eq.(5) 



C L1 = C H1 = 0, C L2 = ^4-e^, C H2 = \f^ 2v , (7) 
where /g is the Fourier transform of wave packet (6) 



5 L = d 2 + ih(-ji + 2 l2 )t/m, 
S H = d 2 - ih(-yi + 2j 2 )t/m, 
(3 = (ad + ix) 2 -y 2 . 



(12) 



Two terms in Eq.(10), (11) labelled by L and H in- 
dices show the contribution of light and heavy holes cor- 
respondingly. In particular the item 1/Sl exp(/3/2(Ji, — 
a 2 /2) in Eq.(10), which can be written as 



Jl = —j= exp 



(k x - k ) 2 d 2 k 2 y d 2 



(8) 



As it follows from Eq.(7), the initial wave packet consists 
of the light-hole and heavy- hole states, and the weight 
of the light hole states is three time greater than that of 
heavy holes. 

Substituting Eqs.(7), (8) into Eq.(5) and performing 
integration we finally have the components of wave func- 
tion *(r,t) = (■9 1 (r,t),^ 2 (r,t),^ 3 (r,t),^ 4 (r,t)) T : 



* 2 (r,i) = * 4 (r,i) = 0, 



de~ a2 ' 2 
4^ 



(9) 



(10) 



V3de~ a2 / 2 a 2 d 2 + (x + iyf 
3(r ' t)_ 4^ (ad + y) 2 + x 2 X 

0/2Sh (1 _ _&Z.\ 
1 (3 } 



_ (1 _2& } _e_ 
Sl K ' S H 



(11) 



where 



exp 



exp(/?/2<5 L - a 2 jl) 
Sl 

jx 2 + y 2 ) + 2iad 3 (x + ggj-71 + 2 72 )) 
2(d 2 +^(-71 + 2 72 )i/m) 

x(d 2 + ^(- 7l + 27 2 )t/m)" 1 , 



reproduces the time evolution of Gaussian wave packet 
for spinless particle with energy El = + 2"f 2 )k 2 . 

In Fig.l and Fig. 2 we represent the hole density for 
initial wave packet, Eq.(6) for 71 = 7.65 and 72 = 2.41, 
which correspond to GaAs, with width d = 3 • 10 _6 cm 
and ko = 2 • 10 6 cm _1 at the moment t — 0.2 (in units of 



t = md 2 /2~f 2 K) (Fig.l); with width d = 3 • 10" 



and 



feo = at the moment t — A (in units of to) (Fig. 2). 

It is not difficult to see in Fig.l and Fig. 2 that the pa- 
rameter a — k d regulates the character of wave packet 
evolution. For the case a >> 1 (see Fig.l) the initial 
wave packets split into two parts and their evolution is 
accompanied by the Zitterbewegung phenomenon. This 
dynamics is a result of interference of the hole states lay- 
ing near the point ko in momentum space. The splitting 
of the wave packets into two parts which propagate with 
unequal group velocity and have different angular mo- 
mentum polarizations, appears due to the presence of 
the light and heavy holes states in the expansion the ini- 
tial wave packet. As usual, the splitting is accompanied 



4 




FIG. 2: (Color online) The full hole density at the moment 
t — 4 (in the units of to) for 71 = 7.65 and 72 = 2.41 for initial 
wave packet, Eq.(6) with width d = 3 ■ W~ 6 cm and ko = 0. 

by the packet broadening which appears due to the effect 
of dispersion. 

The ZB is a result of interference and a complex space- 
time evolution of two parts of wave packet. In another 
extreme case, when the inequality a << 1 takes place, 
the oscillations of the probability density at the circum- 
ference of the wave packet are closely related to the inter- 
ference of heavy- and light-hole states located in vicinity 
of the point k = in the momentum space (see Fig. 2). 
This interference gives rise to the formation of the circu- 
lar ripples around the packet center. The local period of 
these oscillations of probability density is not constant, 
it decreases with increasing radius of the ripple and time. 
This conclusion can be illustrated by the analytical ex- 
pression for the density probabilities for a = 0. For ex- 
ample, the analytical expression for the first component 
of wave function is 




at = 15 
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a = 0.1 



2 4 5 8 10 

FIG. 3: (Color online) The average coordinate y(t) of the 
packet center for different cases of value a = kod. 

representation can be obtained from Eqs.(4), (5), (6). 
After that the usual definition 

4 f d 
Xi(t) = J2 dkC*{k,t)i—C a (k,t), i = 1,2 (13) 

a— 1 1 

leads to 

x{t) = — (72-7i)«, (14) 



1 



9 exp(- 



d 2 (r+w£t 2 ) 



Tffi) ex P(" 



1 + u 2 t 2 



+6 



d 2 (l + LJ 2 H t 2 ) 

1 + uJ 2 H t 2 



(1 + UJ L LU H t 2 )x 



16nd 2 1 

( r 2 r 2 \ 

ex Pl~2rf2( 1+a) ^ 4 2) - 2d 2 (l+aj' A H t 2 )' 

1 + U) 2 L L0 2 H t A + LU 2 t 2 + LU 2 H t 2 

x cos(£(f, r)) + (uj l - u H )t sin(£(t, r)) 



where lu l = ft(— 71 + 2^ 2 )/d 2 m, uj h = —ft (71 +2^ 2 )/d 2 m, 
n = TT^ + TT^' and €(*,r) = ^n. This expres- 
sion contains oscillating terms which are proportional to 
functions cos(£(i,r)) and sm(£(i,r)), describing the rip- 
ples structure at Fig. 2. It should be noted, however, that 
these oscillations do not provoke the ZB of average coor- 
dinate of packet center. 

To analyze the motion of the packet center we have to 
find the average value of the position operator. To do 
it, let us use the momentum representation. The com- 
ponents C a (k,t) (a = 1, ...4) of wave function (5) in this 
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4k~ 



exp( 



{at /to) 



1 + (t/t y 



) cos( 



a 2 t/t 
1 + (t/t ) 



r)-l 



(15) 

We see that the center of wave packet moves along x 
direction with constant velocity Vq x — ^(72— 7i)- This 
expression also can be found from the equation 



Va. 



dEr 



Wl + -V H , 



(16) 



where V^h) — j~^jr^'\k 3; =k are the velocities of light 
(heavy) holes along x axis at k x = ko and E L ^ H ^ are 
determined by Eq.(2). The coefficients 3/4 and 1/4 cor- 
respond to relative parts of L— and H— holes. The mo- 
tion of the wave packet center along x is accompanied 
by the oscillations of the packet center at a > 1 in a 
perpendicular direction or Zitterbewegung (Eq.(15)). It 
is easy to see that such trembling motion has a transient 
character as for other systems with two-band structures 
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(Fig.3 for a = 15). As it follows from Eqs.(14), (15) for 
a given initial polarization of the wave packet the ZB 
occurs in the direction perpendicular to the initial mo- 
mentum pq x — hko- At large enough time t » to the 
packet center shifts in y direction at the values of 



yo 



4fc 



(e" a -1) 



where Sl, 5h and are determined by Eq.(12). The 
results of the calculation of average values of x and y for 
this polarization are 



x(t) = V 0x t + x Z B(t), (21) 
where the constant velocity is defined by the expression 



in accordance with Eq.(15). The Fig.3 for a — 0.1 obvi- 
ously demonstrates the monotonic behavior of y(t). 

ii) As a second example we compute the free space-time 
evolution of the Gaussian packet 



4-(f,0) 



exp(- 



2iP 



(17) 



This expression is obtained as a superposition of two 
components, corresponding to the light and heavy holes. 
It should be mentioned that the relative weight of the 
light and heavy holes is not constant in momentum space 
in contrast to the previous case, Eq.(6). Really, as it 
follows from Eqs.(5), (17) the amplitudes Or,(#Y, are 



C m = Chi = 0, C L2 = l-foyfap* + 1), 
^/2 



C H 2 = ^fh(e^-V3). 



(18) 



Substituting Eq.(18) into Eq.(5) we obtain the compo- 
nents of wave function 



*2 (»*,*) = *4(r ? ,i) = 0, 



T ^ , V2 d a 2 /2 f3e^ 2 ^ 



P f3/25 H 



V3de~ a2 / 2 a 2 d 2 + (x- iyf ( e^ 2 ^ 
4%/2^ {ad - y) 2 + x 2 y 6 L 

2S L ^ 2S H \ 

' s H [ p >y 



, _i . V2 d 2 /2 /e^ 3e^ 25 " 

# (f t) = =e _a /2 1 ; 

31 ' ' 8 V 5l Sh 

s/Zde~ a2 / 2 a 2 d 2 + {x + iy) 2 ( e^ 2 ^ 
4V2^ (ad + y) 2 + x 2 { S L 



(19) 



(1 



25 L 2S H \ 
' ' S H 1 >V 



(20) 



V 0x = 



hko 



r [ 1 1 e" a 

■7i+ 72V 3 1 ? + t 



, (22) 



and XzB{t) describes the oscillatory motion (Zittebewe- 
gung) 



xzb it) 



V3 



exp(- 



aH2/t2 ° ) ~in( aH/tl 
1 + H l + t 2 /t 2 J to 




o , * 

X 



l + t 2 /t 2 ' 2k a 



V3 
4fc 



-(1- 



I exp(- 



a 2 t 2 /t 2 



cos(- 



(23) 



a 2 t/t 



l + t 2 /t 2 / y l + t 2 /t 2 



f 



a 2 t 



a 2 t 2 /t 2 . a 2 t/t 

exp( "TTW^ )sm( T+l7^ 



(24) 



The space-time diagram of this solution, shown in Fig. 4 
again shows the Zitterbewegung of the position's mean 
values. It is interesting to stress that for the packet po- 
larization determined by Eq.(17) the oscillatory behav- 
ior occurs in both x and y directions. The amplitude 
of oscillations in x direction is much less than that in y 
direction. 

As it follows from Eqs.(15), (23), (24) the character 
of trembling motion strongly depends on the parameter 
a = kod as for the other system^iL Really, it is easy to 
see that oscillations of r(t) (which are connected with the 



terms cos( 



art/to 



)) occur if a 2 > it. How- 



ever the value a should not be too large since the ampli- 
tude of oscillations (~ cxp(— ,,J,A )) decreases rapidly 

l + t /tg 

when a increase. 

Notice also that if we consider the plane wave that 
corresponds to d — > oo (i.e. a — > oo and to — > oo) then 
we obtain from Eqs.(14), (15) 



x(t) = — ( 72 - 7i)i, 
m 

3 / r 272^fcgt . 
4fco \ m 



(25) 



and correspondingly for the second polarization (Eq.(17)) 
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FIG. 4: (Color online) The space-time diagram of the packet 
center for the initial wave packet, Eq. (17) with a — 15. 



:(*) = — (V*ys -<*)*, 
to 

— 1 - cos( 

4fcfJ \ TO / 



(26) 



iii) Let us consider now the third example when the 
initial wave packet has the following form: 



*(r,0) 



cxp( 



2(1- 



+ ik x) 



(27) 



Using expression (5) and specifying the coefficients 
^L(H),i{k), we find the components of the wave function 
at arbitrary time: 



(31) 



It is easy to see that for this polarization the time 
dependence of the full probability density is similar to 
the evolution in the case i) in sprite of that here all four 
components of wave function are not equal to zero. Here 
we note only that the average coordinate x(t) — — hk ^ ri t, 
and y(t) has the same form as in the case i), Eq.(15). 
But as we will see below the spin density and average 
spin evolution in this case differ qualitatively from those 
in cases i) and ii). 



III. ALTERNATING SPIN - 3/2 
POLARIZATION 

As was shown firstly in Ref. [19], the spin dynamics of 
spin - 3/2 hole is qualitatively different from those of spin 

- 1/2 electron systems. In contrast to spin - 1/2 electron 
systems, neither the magnitude nor the orientation of 
the angular momentum S are conserved. In spite of the 
fact that the coupling of spin and orbital degrees of free- 
dom cannot be written in terms of external or effective 
magnetic field, the hole spin polarization, as well as the 
higher-order multipoles demonstrate nontrivial dynamics 
in time and can precess due to the interference of heavy- 
and light- hole states^. In the Ref. [19] the main re- 
sults concerning the alternating spin polarization in spin 

- 3/2 hole systems were obtained only for the initial states 
represented by the plane waves. Below we describe and 
visualize the spin dynamics of hole wave packets with dif- 
ferent initial polarizations. For this purpose we consider 
the time evolution of components of the effective spin 
density Si(r,i) and their average values Si(t) 



(32) 



(33) 



*l(T*t) 



#2(r,t) = i 



de~ a ^ 2 



4V27T 



de- fl2 / 2 



4V27T 



e /3/2S L 3 e /3/2<5fr 



4\/27r (a 2 d + y) 2 + x 2 

eP/ 2S - _ 26l e?/ 2S » _ 25^, 
Sl [ P ' S H [ (3 ' 



where *f/(r,t) is a four-component wave function and 

(28) \Er+(r, t) denotes the Hermitian conjugated wave func- 
tion. 

Firstly, we consider the wave packet, corresponding to 
the hole spin oriented initially along z direction: S z = |, 

(29) Eq.(6). In this case it is easy to see that spin den- 
sities are equal to S x (r,t) = S y (r,t) = 0, S z (r,t) = 
||^I'i(r, t)\ 2 — i|5 , 3(r, i)| 2 , where the components of wave 
function ^i(f,t) and \&3(r,i) are determined by Eqs.(9)- 
(11). 

To illustrate the evolution of the hole spin density 
S z (f,t) we plot this function in Fig. 5 for the moment 

(30) °f t ne time t = 15 (in units of i ) and for wave packet 



parameters d = 10 cm and ko 



10 5 



Just as 



FIG. 5: (Color online) The angular momentum density 
S z (r,t) for initial wave packet, Eq.(6) with ko = 10 5 cm -1 , 
d — 10 -6 cm, a = 0.1 at t=15 (in units of to)- 



for hole probability density the space-time evolution of 
the spin density strongly depends on the parameter a. 
In particular, for small values of a as was demonstrated 
in Sec. II the splitting of wave packet is absent and os- 
cillations of the probability density (ripples) arises on its 
periphery. Similar behavior of spin density S z (r, t) for 
a = 0.1 one can observe in Fig. 5. 

Using Eq. (33) it is not difficult to find the average value 
S z (t) for given initial polarization 



FIG. 6: (Color online) The average value of angular momen- 
tum Sz(t) for initial spinor, Eq.(6) for a = 6. 



the components of the wave function at arbitrary time 
Eqs.(28), (29), (30), (31) one can obtain the hole _spin 

density S(r, t) and the components of average spin S(t) 



s x (t) 



V3 a/3 
x cos( 



a 2 {t/t a ) 2 
1 - exp(- , ' \ 9 )x 



1 + (t/t ) 2 



a 2 t/ta 



1 + (t/toY 



(36) 



cos( 



a 2 t/t 

i + (t/t y- 



( a 2 (t /toy s 

i + (t/t ) 2 

t . aH/tp . 

■ Sin ; — ; 777 



(34) 



This expression obviously demonstrates that S z (t) ex- 
periences the transient oscillations for the parameter 
a >> 1 (see Fig. 6). For t/to >> 1 S z (t) streams to the 
constant values 3/4. Note that the length of spin vector 
in this case is not constant. As it follows from Eq.(34) in 
the limiting case d — > oo, i.e. for the plane wave, S z is 



s z {t) 



2 l2 fik 2t 



(35) 



This result coincides with that in Ref. [19]. 

For the second example presented in Sec. II, Eq.(17), 
the initial values of effective spin are 5^(0) = 5^(0), 
<Sz(0) = 1/2, and space-time evolution of spin density 
components Si(r,t) are similar to that discussed above. 

Note that for these two cases the initial_ vector 5(0) is 
perpendicular to the average momentum p = (hko, 0, 0). 
But as was revealed in Ref. [19], the most interesting spin 
dynamics characterized by hole spin precession occurs 
when p(0) is not perpendicular or parallel to p . 

Such initial condition is realized in particular for the 
wave packet given by Eq.(27). Using the expressions for 



Sy(t) = 
t 



V3 ex P("T 



g (t/*o) 



1 + (t/t 
a 2 t/t 



cos( 

t y l + (t/t ) 2 



x sin( 



s z (t) = - 

aH/to 




sm( TTW )+ 

a 2 (t/t ) 2 



4a 2 ^ v 1 + it/to) 2 ' 



1 + [t/to) 



3ex P (-^W^ 



(37) 



1 + (t/h) 



x I cos d , uu yd --Mtttttt™) 



1 + (t/to) 



t 

to 



a 2 t/t 



I + (t/to) 2 



(38) 



These expressions become more simple in the limit case 
d — > oo, i.e. for the plain wave propagating along the x 
axis with wave vector ko 



S x (t) = — — , S y (t)= — —sin — — — , 



s z(t) = - \ 1 + 3cos(— — 2_) 



(39) 
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FIG. 7: (Color online) The transient precession of average an- 
gular moment vector S(t) (in plane (S y ,S z )) about the vec- 
tor fco || Ox for wave packet described by Eq.(27) (a — 10). 
Dashed line in accordance with Eq.(40) corresponds to the 
hole spin precession for the plane wave. 

Last formulas clearly demonstrate the precession of 
vector S(t) about the wave vector fc : the end of vec- 
tor S(t) describes ellipsis on the plane (S y , S z ) 

<5#£ + ^-y 4 > 2 =l, (40) 

where m = v3/2, n = 3/4. 

In our case for the wave packet of finite width the time 
evolution of the effective spin becomes more complicated 
as it follows from Eqs.(36)-(38). Firstly, the component 
S x (t) depends on time, although this dependence is weak 
enough both for a » 1 and a « 1. Secondly, the "pre- 
cession" in (S y , S z ) plane acquires the transient character 
which is demonstrated in Fig. 7. As one can see from Fig. 7 
the component S y (0) = S y (oo) — and the component 
S z changes its values from S z (0) — 1 to S z (oo) — 1/4. 

To illustrate the spin density we plot the components 
S z (f, t) and S y (f, t) in Fig. 8 for the moments of the time 
t = 0.1 and t = 0.5 (in units of to) correspondingly and 
for a = 6. The Fig. 8(b) demonstrates clearly the split- 
ting of the angular momentum density into two parts sim- 
ilar to probability density discussed in Sec. II for a >> 1. 
The angular momentum density changes its sign within 
each part that looks like as multipole structure. 



IV. CONCLUSION 

Using the Luttinger isotropic model we have studied 
the dynamics of wave packets in hole 3/2 effective spin 
systems. The analytical expressions for the four com- 
ponent wave function are obtained and the probability 
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FIG. 8: (Color online) The angular momentum density 
S z {r,t) and S y (r,t) for initial wave packet, Eq.(27) for fco = 
2- 10 6 cm" 1 , d = 3- 10" 6 cm, a = 6 at t = 0.1 (a) and t = 0.5 
(b) (in units of to). 

densities as well as the effective spin densities are visual- 
ized for different initial spin polarizations. Two different 
types of the spatial evolution of wave packets are found. 
In the case when the parameter a » 1 the initially lo- 
calized packet splits into two parts propagating with dif- 
ferent group velocities and having various effective spin 
polarizations. Similarly to the wave packet evolution in 
systems with spin-orbit coupling^ this evolution is ac- 
companied by the oscillations of packet center or Zitter- 
bewegung. These oscillations quickly fade away when the 
distance between two parts of split packet exceeds its ini- 
tial value. Another scenario of the evolution is realized 
for a small value of parameter a, namely when a « 1. 
In this case the probability density remains almost ax- 
ial symmetric with time, but the ripples appear at the 
circumference of the density distribution. These ripples 
arise owing to the interference of heavy- and light-hole 
states located in the vicinity of the point k — in mo- 
mentum space. It should be noted that this interference 
does not stimulate the oscillations of packet center f(t) 
or ZB. 

It is known that the spin dynamics of the holes, de- 
scribed by an effective spin - 3/2, differs qualitatively 
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from that for spin - 1/2 electrons. The main and striking 
peculiarity of hole spin dynamics is the nontrivial peri- 
odic motion of average spin vector in spin space, which 
can be viewed as a precession. We have studied in detail 
such unusual behavior for the hole wave packets. We have 
shown that finite width of wave packet is responsible for 
the transient character of spin precession. The analyti- 
cal and numerical study of the influence of a parameter 
on the angular momentum density and average spin vec- 
tor is presented. For the case when a » 1 the density 
of angular momentum splits into two parts having multi- 
pole polarization, but for a << 1 the angular momentum 
density preserves almost the initial cylindrical symmetry. 

The experimental observation of trembling motion of 
f(t) and spin precession in crystalline solids is a long 
standing problem of condensed matter physics. The pos- 
sible methods of producing of electron and hole wave 
packets as well as the observation of ZB effects was an- 
alyzed in numerous of recent publications. As was dis- 
cussed in Rcf. [20], the hole packets can be excited opti- 
cally in nondegenerate 3D hole gas by using a polarized 
laser beam. If laser beam has a finite size the optically ex- 
ited wave packet will consist of states of light and heavy 
holes lying in the range of dr 1 in vicinity of the point 
k = and can have a definite spin polarization. To ob- 
tain the wishful average momentum ko the electric field 
pulse with duration of some picosecond can be applied 
after excitation. During the process of optical excita- 
tion electrons will excited also, but the electric field will 
separate positive and negative charges. The ZB evolu- 
tion to be robust against disorder provided ujzb > 1 /Ts 



and ujzb > l/ r p where r s and t p are spin and momen- 
tum relaxation times correspondingly. As was verified 
experimentally (e.g. in GaAs and InAs ts is of the or- 
der 10 — 100 ps) these conditions can be fulfilled in clear 
semiconductors . 

An interesting method of observation of trembling mo- 
tion of charged particles in solids was proposed by Rusin 
and Zavadski in Ref. [15]. They note that the typical 
period of ZB oscillation is comparable with the period 
of Bloch oscillations, and so the method which early was 
used for observation of Bloch oscillation in superlatticies 
at the presence of electric field^ can be used also for the 
ZB observation. The results of Sec. II and Sec. Ill of 
this work show that this method can be used also for ZB 
observation in Luttinger spin - 3/2 system. 

The information about the evolution of hole wave pack- 
ets is important in the context of application for analysis 
of functioning of different spintronic devices (e.g. Datta 
- Das hole spin effect transistor—) and for understanding 
the transport phenomena and spin related dynamics. 
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